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We calculate the density of stationary points and minima of a A*' ^ 1 dimensional Gaussian energy land- 
scape. We use it to show that the point of zero-temperature replica symmetry breaking in the equilibrium 
statistical mechanics of a particle placed in such a landscape in a spherical box of size L — Ry/N corresponds 
to the onset of exponential in A'^ growth of the cumulative number of stationary points, but not necessarily 
the minima. For finite temperatures we construct a simple variational upper bound on the true free energy 
of the R — <x version of the problem and show that this approximation is able to recover the position of the 
whole de-Almeida-Thouless line. 

■PACS: 05.40.-a, 75.10.Nr 



Statistical properties of stationary points (extrema 
,as well as saddles) of random high-dimensional land- 
'scapes continue to attract considerable interest in the 
.physical and mathematical literature, both analytically 
'and numerically. Motivations come from fields of re- 
isearch as diverse as the study of disordered and glassy 
systems and string theory - see [2 [H [3] and references 
'therein. In particular, considerable effort [1] has gone 
into trying to understand what changes the intricate sce- 
nario of ergodicity breaking typical for mean-field spin 
iglass models implies in the statistics of the associated 
free energy landscape (the so-called TAP 5 variational 
functional). 

One of the simplest models known to display many 
non-trivial hallmarks of glassy dynamics and thermody- 
namics is the so-called toy model of a classical particle 
in a box of characteristic size L filled with a random 
potential. The position of the particle is characterized 
by the coordinate vector x = (xi, ...,xn), |x| < L, and 
the standard choice of the energy function is 



k=l 



N 



(1) 



with /i > 0. The random Gaussian-distributed potential 
y(x) is assumed to have zero mean, and its covariance 
is chosen to have a form ensuring stationarity, isotropy, 
and a well-defined large A^— limit: 



(F(xi) F(x2))^ = iV/( 2^(xi-x2) 



(2) 



Similarly, the natural scaling of the sample size is 
L = R^/N 6\. Here and henceforth the notation {■ ■ ■)y 
stands for averaging over the random potential V{x.). 

In the limit when both N oo and i? — > cxd the 
thermodynamics of the model has enjoyed a long his- 
tory of research starting from early works [7], [8], and [9]. 
A detailed discussion of the limit A^ — > oo for the fixed 
box size R < oo can be found in a recent paper [6]. The 
paper [TT] investigates the related problem for finite di- 
mensions A^ < oo in a broad physical context. 

The calculation in [IJ revealed that the replica sym- 
metry breaking (interpreted in the standard way as er- 
godicity breaking) in the zero-temperature limit of the 
R = oo version of the model is accompanied by the 
emergence of an exponentially large total number of sta- 
tionary points in the energy landscape. As naively one 
may expect that the low-temperature thermodynamics 
should be dominated by minima rather than by the to- 
tality of stationary points, the issue deserves further 
investigation, in particular in the general R < oo case. 

The goal of the present paper is to compute the den- 
sity of all stationary points as well as of only minima 
at a given value of the potential, hence energy, in the 
limit A^ ^ 1 (in the simplest case A^ = 1 similar ques- 
tions were addressed in [IHI)- Our analysis reveals that 
generally the domain of existence of the glassy phase 
with broken ergodicity (at zero temperature T) turns 
out to be associated with the existence of exponentially 
many stationary points in the energy landscape, but not 
necessarily exponentially many minima. In an attempt 
to extend our considerations to finite temperatures we 
also construct a simple variational functional providing 



1 



2 



an upper bound on the true free energy of the R = oo 
version of the problem. Surprisingly, counting station- 
ary points in this simple-minded approximation is able 
to capture such a nontrivial feature as the precise posi- 
tion of the de-Almeida-Thouless line in the whole (/i, T) 
plane. 

Our aim is to calculate the mean probability den- 
sity ps{V,x) of the value V of the potential V{x) taken 
over all stationary points of the landscape 7i(x) around 
the position x, with the stationary points of any index 
counted with equal weight. We will also be able to ad- 
dress its counterpart, p,n{V,x.), where only minima of 
7i(x) are taken into account. It is convenient to define 
the joint probability density !F{V, x, K) of the scalar ar- 
gument the iV— component vector argument x and 
the N X N real symmetric matrix argument K as 



d 



(3) 



in terms of which the density Ps(V, x) can be expressed 
as [J 

(F, x) = j \ det (^i/jv +K)\T{V,^,k) dk (4) 

where /jv stands for N x N identity matrix. The den- 
sity of minima pm {V, x) is given by a similar integral, 
but with the integrand containing the (matrix) step- 
function factor 6 (^filN + ensuring that all the eigen- 
values of the matrix pl^ + K are positive. 

Introducing the Fourier integral representation for 
each of the delta-functional measures (of scalar, vector, 
or matrix argument) in Eq.Q facilitates performing the 
ensemble average explicitly in view of the Gaussian na- 
ture of the potential Eq.(IT]) and its stationarity, Eq.Q. 
After some straightforward but long manipulations we 
find: 



J-"(y, X, K) cx exp 



2/'(0) 2iV/(0)J 



N +2 N 



:Tr[i5--(9*+^f)/«] 



dt (5) 



Here, and henceforth, we systematically disregard var- 
ious multiplicative constant factors for the sake of 
brevity and use the following notations: 



Pc 



5^ = /"(0)-^>0. 



(6) 



Substituting Eq.(l5]) back into Eq.(|3]), and introducing 
the notations 

(0) 



we find that the density function factorizes as 

p,(y,x) -AA,po(x)p.(y), (7) 

where 



po(x) = e^FIoT, ps{V) = e ^w/co) 



N J.2 

e'^* Vs{t,V)dt 



(8) 



with 



{t, V) = f\ det (^p,f{t, V)iN + i?) |e 



dH. 



(9) 

and the factor Ms accounts for the overall normalisation. 
Analogously, the density of minima Pm(V, x) is given by 
a very similar expression, the only difference being that 
Eq.® is replaced by 



X det (^pefit,V) 



N 



H]9 



In+H] dH. 



(10) 



Our main interest is to extract the leading asymp- 
totics of these expressions in the limit ^ 1. To this 
end, for the analysis of Ps{V) we use the results from 
[l], where it was demonstrated that essentially Ds{t, V) 
can be replaced in the large N limit with 

I?,(i,l^)cxexpAr$(s), p,f{t,V)/pcr (11) 

where the function ^(s) ~ ^{—s) is given explicitly by 

„2 



n-s > 0) = --9(3-2) 



In 



(12)- 

Perhaps, the shortest way to understand the above 
asymptotics is to notice that actually, 

$W= J \n\s + X\i,sc{\)dX, :/,,(A) = — 

' . . ^ ^^^^ 

where the integral is understood in the sense of the prin- 
cipal value. The above formula can be verified by ap- 
plying ideas from statistical mechanics to the evaluation 
of Ds{t,V). A rigorous exposition of the method and 
further references can be found in [T2j, and a related 
informal treatment is given in a recent insightful paper 

m- 

To ensure a well-defined large A^ limit we rescale 
V — > NV, and evaluating the integral over t with the 
Laplace method we find that Ps{V) cx exp [— A^5(V^)], 
with 



H = K- 



^ f'(0) V\ f 



/(O) 



V 



mmf 



|-$(^^^(t, V)/p,r) 



(14) 



The density of stationary points 
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Correspondingly, in the large N limit the density func- 
tion PsiV) displays a pronounced maximum around the 
most probable value V = = argminiS(V^), which can 
be found from the system of equations: 



/'(O) d<i>{.3) . 



9 d<i>{s). 



(15) 



j-l'C.f GiS f-^cT" CIS 

where we have used the notation 

= [/^+.9^* + 7^ ^*]/a*c. ^ [/.+-^ K]/Mc. , (16) 

the second equality following from the obvious relation 
t* = 5K//'(0) implied by Eg. p^ . and the definition 
Eq.®. Moreover, for > Eq.dlSl) imphes that 



ds 



(17) 



Assuming first < < 2, the first of the equa- 
tions pS]) together with Eas. (|16ll7p immediately yields 
y* = ^f'{0)/fM^^ and s* = 2^//icr, the solution being 
compatible with our assumption for < ^ < ficr- On 
the other hand, assuming > 2 the first of the equa- 
tions (fT^l) solved together with Eq. pT)) yields the rela- 
tion 

-1 



s ^ 



/'(o) 



y* 



/'(O) 



v; 



(18) 



which is consistent with our assumption. Substituting 
for from Eg. p^ one immediately finds Vt = f'{0)/fi. 
Summarizing, we see that the most probable value of 
the potential ^(x)) at a stationary point of the energy 
surface 7i(x) is given by 



0ifl - /i,,) + 0(/icr - /i) (19) 



The corresponding values of S{V^,) can be found from 
Eq.lll]): 



5(K 



2 



In 



(20) 



ordinate integration is easily performed by the Laplace 
method yielding 



\yi\<RVN 



/9o(x) dx a 



2 I In (27ri?2) + 1 - |^ , R<R 



(21) 



Collecting all exponential factors from Eqs. ()20l2ip and 
also restoring the overall normalisation from the require- 
ment lim^^oo(T^s)y = 1, one finds that the associated 
cumulative complexity of stationary points Y,s{p-,R) = 
limjv— oo jq In is given by 



\9{R^ - R) 



Rl 



1-^ 
Rl 



< (22) 



as long as fi > ^cr, whereas for /i < ficr the complexity 
Ss(/^, i?) is given by 



- 2 

y - i 



-H 
ln(f 



for 

R < Rn 



R> R, 



where introduced an important lengthscale, cf. [6] 

Rcr = vuwnw) ■ 



(24) 



The calculation of the density of minima PmiV) pro- 
ceeds in a similar way. The required large N asymp- 
totics of Vjnit, V) in Eg. pH)) can be extracted from the 
results of a recent paper ^13j. and is given by 

Vm{t, V) cx exp ( -— ^(2 " s)g{s) + iV0(s) j , (25) 

where s has the same meaning as in Eg. pTj) . 0(s) is 
given by the right-hand side of Eq. p^ for s > 2, and 



- 216 



72s2 _ s4 



(s3 + 30s)Vl2 + s2 
' s + Vs^ + 12\ 



log 



6 



(26) 



The mean total number (#s)y of stationary points 
is obviously obtained from the density Ps{V, x) in Eq.Q 
by integrating it across the range of coordinates inside 
the sample |x| < i = R\/N and over all possible values 
of the potential V . The integral over V is dominated by 
the most probable value found above, whereas the co- 



When calculating Pm{V) we again rescale V —> NV , 
subdivide the integration range into two pieces s € 
(—00,2] and s G [2, 00) and evaluate two resulting inte- 
grals by the Laplace method in the large N limit. The 
first integral Ii {V) is dominated by the point where Q (s) 
is minimized, i.e. by the vicinity of s = 2, the endpoint 



4 



of the integration range, where actuaUy Q = 0. Remem- 
bering that 4>{2) — 1, we find after some manipulations 
Ii{V) cx exp-NSi{V), where 



1 



25V(0) 



JL 

f-^cr 



1. (27) 



Evidently, /i(V) is strongly peaked about the value Vi, 
where 

(28) 

In the second integral, hiV), the arguments lead- 
ing to Eas. (|19l20p show that for any n it is dominated 
by the same value V = V^, = f {0} / fj. as in the /i > /icr 
case in (fT9l) . and is given by hiV) cx exp —NS{V^,) with 
iS(T4) = —1/2 -|- In {^cr/f^) (see (|20)) ). It is easy to ver- 
ify that iSi(Vi) > iS(t4) for > ^cr, and for /i < ficr 
the inequality is reversed, giving the leading large N 
asymptotics as 



Pm{V) cx 



exp— A^5(V*), fi > ^cr 
exp-NSi{Vi), fj.< jJc 



(29) 



We thus find that the associated complexity of the 
minima, given by ^mi^J■, R) = limAr^oo jj In (#m)v', co- 
incides with Ss(M7 R) in the range fi > ficr- For fi < fj,cr 
the counterpart of ([23]) is however given by a different 
expression: 



1 



> (30) 



as long as R> Rf^, whereas for i? < i?^ it is given by 



temperature the phase with broken ergodicity covers 
the strip Rcr < R < oo, < /i < in the {R, fJ-) 
plane. As is easy to infer from the last line in ex- 
pression Eq. (|^5)) . the above strip is exactly the domain 
supporting positive cumulative complexity > 0. At 
the same time, analysis of (|30|) shows that Em > for 
any < fi < ficr only when the sample radius sat- 
isfies R > Rm — eRcr > Rcr- On the other hand, 
if Rcr < R < Rm then for the values of fj, satisfying 
< /i < /i_ < j-icr with 



2fic 



1 + R^/Rlr 



1 



1 - 



1 + R^/Rlr R. 



2 R 

(32) 

the complexity of minima S„j is negative, although we 
have seen that > in this range of parameters. 
We thus conclude that the replica symmetry - hence, 
ergodicity,- breaking implies that the landscape sup- 
ports exponentially many stationary points, but not 
necessarily minima. 

It is natural to try to understand whether similar 
relations exist at finite temperatures for the free energy 
landscape. For the present model the full TAP varia- 
tional functional is as yet unknown, and constructing 
it we consider as an interesting problem for a future 
work. To this end, we propose below a simple varia- 
tional functional which provides only an upper bound 
for the true free energy of our problem. Surprisingly 
enough, we find that such a functional already contains 
enough information to reproduce the position of the de- 
Almeida-Thouless line in the whole (/x, T) plane. 

Concentrating, for simplicity, on the infinite sam- 
ple R — oo case, we employ for our goal the Gibbs- 
Bogoliubov-FeynmanHJ inequality: 



F{H) < F{Ha) + {H-H, 



(33) 



2- 



R^jj^ 
Rlf^l 



In 



R^ 



- 2- 



(31) 



Let us first compare ([23)) and ()30l3ip in the infinite- 
sample limit R — oo studied in |Jj. In that case, we see 
that both an exponentially large number of stationary 
points, and exponentially many minima co-exist in pre- 
cisely the same interval < /i < ficr- Moreover, it is 
precisely the same interval of the /i— axis known to cor- 
respond at zero-temperature to the glassy phase with 
broken ergodicity [71 [H [6] . 

Recent analysis of the general R < oo case of our 
model by the replica trick [6J has revealed that at zero 



valid for any two Hamiltonians H and Ha, with F{H) = 
-TlnTre'"/'^ and F{Ha) = -TlnTre"^"/'^ standing 
for the corresponding free energies and 



((■■•))h =Tr (...)e 



/Tre 



(34) 



The role of H in our consideration will be played by 
7i(x) from Eq. ([!]), and we choose Ha in the form 
Ha{x) — §(m — x)^ with the A''— component vector m 
playing the role of a variational parameter. For this 
choice, the inequality Eg. ([55]) takes the form 



F{n) < Fa{m) 



2 T 



m^ + Va{m) (35) 



The density of stationary points 
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where 



K(m) 




(36) 

with dy = Hill "tIt' "^^^ problem of finding the best 
possible approximation to the true free energy in this 
class of trial Hamiltonians obviously amounts to mini- 
mizing Faijna) with respect to the parameter m. Note 
that Vaijoa) is obviously a random Gaussian function 
with zero mean. The corresponding covariance function 
can be calculated using Eq.([2]) and Eg. ([55]) as 



{Va (mi) Va {V^2))v =Nfa\ T^Av^l - m^) 



2N' 



(37) 



where 



cxp<j j 

(38) 

with the function / taken from Eq.([5]). The invariance of 
the integrand in Eq. ensures that the last expression 
indeed depends only on b^. Choosing b = (|6|, 0, . . . , 0), 
and rescaling y (yi, \/N y ) , we can easily evaluate the 
integral over y in the large N limit by Laplace method. 
This immediately leads to a very simple relation: 



hm fa (b^) =f{h^ + - 



(39) 



Now it is evident that the free energy landscape repre- 
sented by the variational functional Fa{m) in Eg. ([55)1 
is very similar to that of the zero-temperature energy 
landscape 7i(x) from Eq.(IT]). In particular, the func- 
tional Fa{ni) should possess exponentially many sta- 
tionary points, and also minima, as long as 



(40) 



where we recall the definition Eq.®. The condition 
Ea. l|40p coincides precisely with one describing the posi- 
tion of the de-Almeida-Thouless line in the (/i, T) plane 
for the present model [3 [H [6] . 

When the major part of the calculations presented 
in this paper were already completed, we became aware 
of a recent preprint by A J Bray and D S Dean [15] with 
a quite similar scope. The authors restricted their pre- 
sentation to the limiting case = of the same model, 
but were actually able to calculate complexity of critical 
points with a given index. A characteristic lengthscalc 
Lc appearing in their paper is simply related to our 



factor accounting for the difference between their cubic 
and our spherical shape of the sample in the large N 
limit. 
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